RUB-TPII-03/07 



QCD factorization approach for rare B° D* 7 decay 



Nikolai KiveQ 

Institut fur Theoretische Physik II,Ruhr-Universitat Bochum, B-^180 Bochum, Germany 



Abstract 

We present the estimate of the branching ratio for the rare decay B° — » D*j. We use QCD 
factorization approach in order to compute the amplitude of the process. The calculation is carried 
out with the leading order accuracy. The appearing non-perturbative matrix elements have been 
estimated using the large— N c limit and QCD sum rule approach. We obtained that B(B° — » 13*7) — 
1.52 x 10 -7 . Such value of the branching fraction is too small in order to be measured at present 
experiments. 
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1 Introduction 

The different decay processes mediated by b — > cdu quark decay attract a lot of attention of both 
experimentalist and theoreticians. Corresponding hadron decays include various processes like B — > 
D + (it, p, K, ...). There are a lot of experimental results for the different decay modes, see for instance [1] 
and the references there. From the theoretical side the progress in the phcnomcnological description of 
the data can be related with the factorization approach developed last years. The factorization theorems 
for the different decay channels have been discussed in literature [2j [3l |H [5j [6l [7] . In the present paper we 
would like to consider one particular decay mode B — ► D* +7 which remained beyond the considerations 
mentioned above. 

From the experimental point of view the process can be clearly observed due to the higher energy of 
the outgoing photon (E 7 ~ 2.3GeV). The search of this rare decay have already been made by CLEO 8 
and BABAR |9J collaborations. Despite the process has not been observed (B(B° — > D*j) < 2.5 x 10~ D 
[9]) the increasing statistics of the B— factories may provide new opportunities for the better analysis. The 
various existing theoretical models [TU1 [TTl [T2] estimate the branching to be of order of 10~ 6 . Potentially, 
such cross section can be observed despite to the small value and therefore the more qualitative theoretical 
analysis is desirable. 

In present paper we use the factorization technique developed last years for the heavy quarks decays 
in order to derive the leading order factorization formula for the amplitude of the process and estimate 
the branching ratio. Our presentation is organized as follows. Sec. 2 contains the necessary definitions 
and derivation of the leading order factorization calculations. In Sec. 3 we consider the arising soft matrix 
elements and construct the models for these non-perturbative functions using large— N c limit and QCD 
sum rules. This section contains also our main results, the summary and discussions. 

2 The leading order amplitude 

The decay amplitude B°(Pb) — > D* (P[>)j(q) is given by the matrix element 




which is described by the two form factors Fi 2 - Here we accept standard notation a = §- ~ 1/137 
and e* D denotes photon and D— meson polarization vector respectiveljU. The kinematics of the decay 
is very simple. As usual, we choose the frame where B— meson is at rest. Then the the components of 
the momenta read 
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where we introduced the light-cone vectors n, n and for arbitrary vector a one has 
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Figure 1: The leading order diagrams a, b denote the graphs for the form factors a f 
and a nf respectively. The crossed lines denotes the emission of the photon. 



Substituting the numerical values of the heavy meson masses Mo = 2GeV and Mb = 5.28GeV one finds 
Ery ~ 2.3GeV, i.e. the photon energy is quite large. The width is given by 
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Using the experimental constrain for the branching [3] 

B(B° -> ZT°7) < 2.5 x 10~ 5 
and the lifetime t b o — 1.536 x 10~ 12 s one can find for the combination of the form factors in 
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where the coefficient Gp\ / 2ira is introduced for convenience. 

Our task is to compute the form factors F± t 2 i n the limit mb, m c — > oo with m c /mt fixed. The effective 
Hamiltonian in the matrix element ([1]) reads 



H, 



eff 



Gp 

7f 



v cb v: d 



Ci (cb) V - A (du) v ^ A + C 2 (db) v _ A (cu) 



>v-A 



(12) 



where as usually V — A = 7^(1 — 75) and the color indices are not shown explicitly. Let us introduce the 
following parametrization for the amplitude ([T]): 



A D ^ = ^M d MbV2^ G f VcbV* d [ a f + a nf ] , (13) 
where the coefficients a and a are related to the matrix elements of the two operators in the (fl"2"l) : 

a f =i[dxe i M e^( K v'\T[j; m {x) 1 C 2 {db) v ^ A {cu) v _ A )\v) 1 (14) 
a n! =ijdxe t{qx ^ e^(v'\T{j e ; n (x),C 1 (cb) v ^ A (du) v ^ A }\v), (15) 

The meaning of the superscripts "f, nf will be explained below. In these formulas we assume that the 
meson states \v), (v'\ have mass independent HQET normalization. 

In the large mass limit the energy of the photon is also large E 1 ~ itiq — * 00. The emission of such 
higher energy photon is related with short distance subprocess. In some sense, the similar situation is 
encountered in the case of semi-leptonic decay B — * jlu. The difference with respect to our case is in 
the more complicate structure of the matrix element (JXJ) - Consider the simplest diagrams which can 
contribute at the leading order FigfT] The analytical expression for the case a f reads 

a f = AiC 2 Jdxelrt (v 1 \e d c ^P L u d{x) e*k(x, 0) 7ai P l 6(0) + e u d{0)^P L b c(0)^P L A(0, a:)e;«(x)| v) , (16) 

where A(x, y) is the fermion propagator in position space, e Ut d are quark charges, and the hat denotes 
the contractions with the Dirac matrices 0^7^ = a and Pj, = |(1 — 75). In the heavy quark limit one 
performs transition from the QCD heavy quark fields to the HQET fields: 



b(0)~H v , c(0)~h v > 



(17) 
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Hence 



4C 2 z j dx e l{qx) (v' e d h v , >fP L u d(x) £*A(x,Q)^P L H v + e u d>fP L H v h v , T A(0,x)i*u(x) (18) 

e;A(-k + g)ypj f dx (v' \h v , lp P L u [d a (x) (H v ) ] \ v) + 



4iC 2 J dk e d 
AiC 2 I dl e u 



7 "Pl M-l ~ q)i: 



a/3 



dx e" l(fa) («' | [(Ma u^x)} d lp P L H v \v) . (19) 



In the second line we performed the transition to the momentum space, indices a, (3 denote the spinor 
indices, dk = d 4 k/(2ir) 4: . To proceed further we assume that given expression is dominated by the region, 
where the momenta k and I are soft: 



h ~ h ~ A, k 2 ~l A ~ A J 



(20) 



where A ~ Mq — tuq is the soft scale. Then in that region the expressions for the matrix elements 
(i/|...|i>) in (|19p are defined only in terms of the long wave fields and can be understood as soft matrix 
elements. The expressions in the [...] Q| a can be simplified for the large energy (i? 7 ~ tuq): 



tM-k + q)Y l PL = t - l \ k X*) . l^L d S*q^P L 



(— k + q) 2 + ie 



^P L A(-Z ~ <lK = I^Pl 



i(-l-q) 
(-1 - q) 2 +ie 



E* ~ YPl qi 



2(kq) + ie] ' 
—i 

1 mq) + ^e] 



(21) 
(22) 



Substituting these expressions into ([19 
a 1 = AiC 2 J dk 
+ 4iC 2 [ dl 



ie d \elq^P L ] . , 
'-'•■^ ' •"• a ' 3 1 dx e i{hx) (v'\h^P L u [d a (x) H fi ]\v) 



[-2{kq)+ie] 
ie u [rP L qe;] a(3 



m q ) + fe] 



dx e~ 1 ^ («' | [h a u {x)) d^P L H | v) 



(23) 



Performing integrations over dk- and dk± and then over the conjugate variables x+ and x± (and similar 
for the second term with momentum I) we obtain 



C 2 - J dk + { _ 2E ^ k+ + i£] J d(Ai n) lp P L H\v) 



-C 2 - j dl + 



Mi^q £*j"Pr} fd\ 2 



[2E 1 l + + ie] J 2ir 



a + X2 (v' | h lp P L u (A 2 n) d lp P L H\ v) 



(24) 



The formula (|24|) represents the form factor a f as a convolution of the soft light-cone matrix elements 
with the expression which, obviously, is associated with the hard coefficient function. The arguments of 
the fields which are not written explicitly in the eq. ([24|) are set to zero. From the structures of the traces 
one observes that only the combinations antisymmetrical with respect to exchange ji <-> p survive in the 
soft matrix elements. Therefore we define 



S [ d ap] (k+) = AS f ^ e lfc + Al (v 1 \h v , la P L u J(Ai n) lp P L H v \v) , 



(25) 



S^(l+) = AS y ^ e- a ^ (v' | h v , lp P L u (A a n) d 7a P L H v \ v) , 

where symbol "AS" denotes antisymmetrisation with respect to indices {a, p}, for instance 

1 



AS n a ep = - (v 



(26) 



(27) 
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The parametrisation of these functions can be written afH 

Sl ap] (l+) = W f ( l +) AS "° Wh) P ~ ie±p»*%} + AS n a {{e* D ) p + is ±p ^} , (28) 

S l I" ] (k + ) = ±D l (k+) AS n a {(e* D ) p - ie Xp ^} + ^D i {k+) AS n a {(e* D ) p +ie ±pll e% } . (29) 
Then the final result for the form factor reads 

a f = [ {e* ■ e* D ) + ie^^hU^U ie * C * V ' + [ ( e i ' £ h) ~ ^± ^„C 2 U\ (30) 

where we introduced the convolution integrals 

^=r dk+ ^p-, is-r^ZM (3D 

Jo fc + Jo ' + 

The similar calculation for the second form factor a nf provides 

a af = [ ( £ ; • + isf d 2? nf + [ ( £ ; • e* D ) - fcf (ej,),.^),,] i e„C 1 W\ (32) 

where the convolution integrals 

^ = r dk+ ^M , ^=r^ + ^±i (33) 

Jo ^+ Jo '+ 

include the contributions from the different soft matrix elements 

AS / e- u + X2 (v 1 1 h v , lp P L H v d la P L u (A 2 n)| v) 
= i C/ nf (/+) AS n CT {(ej,) p - ie^eS*} + | l> f (/+) AS n CT {(^) p + ie^ej,"} , (34) 

e 4fc+Al («' |^7 ct -Pl^ rf(Ai n) 7p P L u| w) 

= |£> nf (fc+) AS n CT {( £ ^) p - iexp^"} + f £ nf (fc+) AS n CT {(e* D ) p + ie Xplt E%} . (35) 

Let us briefly comment the obtained results. We have performed the calculation only of the leading order 
diagrams. The matrix elements of the non-local four-fermion operators (|25l26p and (|34I35[) consist of the 
product of two field substructures: local one and non-local one. The non-local part is presented by the 
two quark fields separated by light-cone distance. It is clear that such block is not gauge invariant and 
therefore the answer is not complete. To restore the gauge invariance one has to consider the diagrams 
with the emissions of the soft gluons from the active quark. This will restore the gauge link 

E [Ai, A 2 ] = Pexp (*9 J du (Ai - A 2 ) A+ [(uXi + uX 2 ) n]\ (36) 

which connects the fields and therefore completes the definitions of the soft operators. We do not present 
these details because they are standard. One can avoid that using the light-cone gauge A + = 0. Then 
the gauge link (|36|) equals to one and the definitions (|25I26|) and (|34l35p in this case are exact. 

The important question which has to be considered is the exitance of the convolution integrals (|31|) 
and (|33[) . In order to answer it one has to consider the next-to-leading order calculation of the amplitude 
or at least the evolution kernels of the soft operators. Moreover, such calculation is important in order 
to perform the summation of large logarithms which usually appear in the radiative corrections. From 
our calculation we observe that typical virtuality of the hard quark is of order ~ AEj, i.e. we computed 
the leading order contribution to the so-called jet function. The loop corrections contain also corrections 
from the different hard subprocess with the virtualities of order ~ TOq. The presence of the two large 



we use notation ie p ? = \ie pa ^ n p n v . 
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scales unavoidably leads to large logarithms mentioned above. In order to formulate the factorization in 
the general case (i.e. valid to all orders in the QCD perturbation theory) it is convenient to involve the 
technical approach known as soft collinear effective theory (SCET) [T3l I14j . In the present paper we do 
not provide such detailed analysis and restrict our consideration to the phcnomenological estimate of the 
decay width © using the leading order formulas ([50)1 and (|3"2"1) . Below we consider various models for 
the soft matrix elements which we need for the numerical analysis. We shall see that these models are 
in agreement with the factorization, i.e. they have the appropriate end-point behavior which makes the 
convolutions integrals well defined. Of course, this is not a proof but it can be considered as an indication 
that the factorization in the case under consideration is not destroyed by the end-point singularities. 



3 The soft matrix elements and decay width 

Our task is to estimate the non-perturbative matrix elements defined in the previous section. The 
corresponding functions F f ' nl — {?7 f ' nf , D f,nf , U { ' ni , D { ' ni } depend on momentum fraction of the light 
quark velocities v and v', and factorization scale /i^ . In general one can write 

F f,nf = v+F W( k+/v+ , v > +/v+i {v . v ') t MF ) (37) 

The values of the (v ■ v'), v' + are fixed by kinematics and we shall not consider this arguments as an 
arbitrary variables. The factorization scale [If we shall assume to be of order AE 7 ~ 1.5GeV. Usually, 
in that case one has to consider the resummation of the large logarithms which appear in the radiative 
corrections. We do not consider this question in this paper. In future we shall continue to write only one 
argument k + as before to avoid the complexity of the notation. 

Using the time reversal invariance of the strong interactions one can show that the functions F l ' ni are 
real functions. As we shall see later this statement is naturally realized in our models. 

Consider the limit N c — > oo. As one can easily observe 

Ci ~ <D(N°), C 2 ~ OiN- 1 ). (38) 

But for the matrix elements: 

D^U^Nc, D ni ~ U nf ~ A c °. (39) 

Hence both form factors a f,nf are of the same order with respect to large- N c . Note that in our analysis 
it is assumed that we first take the limit itiq — > oo and after that N c — > oo. The conclusion is that 
despite the soft matrix elements have the different order with respect to large- N c we must consider both 
contributions a f and a ni . However the large- N c limit analysis is useful because it allows to estimate the 
contributions to a f . 



3.1 Form factor a f 

The corresponding matrix elements has the factorisable structure and can be approximated at the large- 
N c limit as the product of two matrix elements. We have two non-perturbative functions corresponding 
to non-local d— (|25|) and u— quarks ([26| . For the case of d— quark we can write 

j ^ e l k + Al n p AS (v' \ h v > jx*P L u d (Ai n) lp P L H v \ v) 



^1 e - k+M l(v'\h v , 1±l7 P L u\0 ) (0| d{Xi n) 1+ P L H v \v) 



-A £ *D~)<t F st 



~iF st 4>+{k+, 



'l' F st Mk+h 



(40) 
(41) 



where F st is the static mass-independent decay constant in HQET which is related to the physical constant 
of the heavy meson decay as 



f Q ^M^=F st (l + 0(a s )). 



(42) 



G 



The function <fi + is known as B— meson light-cone distribution amplitude (LCDA ) [15]. Combining ([41]) 
with the parametrization ([29]) one obtains that at the large- N c limit 

D\k + ) = D f (k + ), (43) 
D\K) = 1 -F%<j> + {k + ). (44) 

For the second operator ([26]) one has 

/ "2^ e " li+A2 "" AS ( w '' £ l5l^' 7±^u(A 2 n) J j p P L H v \v) 
a. - x I /■ e -^ + A 2(u / ; ^|^, 7 - Piu( A 2 „) | ) ( | J A M» (45) 

Z7T 



-- F 

2 st 



\{e* D \ 9v{l+)-ie7{e*n)» 9a(1+)) , (46) 



where we introduced the transverse LCD As: 



^ e" 1 l + X2 (v',s* D \h v , 1 x,Plu{M n) |0) = \{e* DA ) a F st /v' + g v (l + /v' + ) 



-~s^{e* D y F st /v' + g A (l+/v' + ). (47) 



These new functions can be related to the LCDA </>+ due to the heavy quark spin- flavor symmetry [171118) . 
The corresponding relation reads, cf.[15 : 

9a(1+) - 9v(l+) = (48) 
Combining ([28]) , J46]) and gg]) one finds 

U { (1 + ) = ^ F^ + (l+/v' + )- (49) 
Hence for the convolution integrals (|3ip we obtain 

, [°° D f {k + ) 1 2 p + (fc + ) F 2 



DC- 



««. / ^^^'^r^Mti.^.?. (5i, 



Substitution of these values into ([30]) gives 



The quantity As is well known from the phenomenology. It was also estimated with the help of sum rules 
[T5], [T6] - For the numerical estimate we accept the value As(lGeV) = 0.35±0.1GeV. For the static decay 
constant we use the value [UlEO] F^lGeV) = 0.35 ± 0.05GeV 3 / 2 and for the coefficient function in the 
effective Hamiltonian ([i"2]) we accept the leading order value Ciimb — 4.8GeV) = —0.268 [21] . Then 

IOV ~ (e* ■ e* D ) (0.97GeV 2 ) + ie ±ap e* a e^ (7.1GeV 2 ) . (53) 
3.2 Form factor a nf 

The two remaining non-perturbative functions D and U which contribute to the a nf can be estimated 
using the method of QCD sum rules. For this purpose consider the following correlation functions (CFs) 

K^(u i u' ) X ) v-v f ) = 1 fdxdye-^+^'y^' (0\T{ r D (y),O°(\),J B (x)}\0), (54) 
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Figure 2: LO diagrams for the perturbative and non-perturbative spectral densities. Diagrams with the 
gluon and quark-gluon operators are not shown. 



where we used following notation 

OaW = \ [Kh1PlH v d{\n) 1+ P L u - h vn+ P L H v d(\n)ilP L u] , (55) 

OuW = \ [KhIPlHv d 1+ P L u(\n) - h vn+ P L H v d^lP L u{\n)] , (56) 

Jd{v) = u(y)j u h v ,(y), Jb{x) = H v (x)ij 5 d(x). (57) 

The index q — u, d is used to specify the non-local structure of the operator. Each CF K± q is parametrized 
by two form factors: 

K™ = (gT -ieT)K q + {gT + ieT) K q (58) 
Saturating the correlation functions with hadron states one obtains for the relevant form factors 

^(i-.KM ^ + "•' (59) 

K * - 4( a-4ta^) ^ + (60) 

where dots denote the contributions from higher resonances and continuum . 

On the other hand for large negative lo, lo' form factors K, K can be computed in Euclidian region. 

K d = ( I ds Pd (s,s',\), K u = ( I ds p u (s,s',X), 

J S — LO J S' — Lj' J S — CO J S' — LO 1 

where spectral densities receive contributions from perturbation theory and from vacuum condensates 

p = p pclt + p coad 

The leading order diagrams for the perturbative and non-perturbative contributions are shown in Fig.l. 
Performing the subtraction of the the continuum contribution (loq is continuum threshold) and introducing 
the Borel transformation with respect to lo and to' one obtains 

l -Fl D{\) e-^ = 4 J ds J ds' e -(^')/2t pd{s , s >, A) , 

l -FlU{\) e-^ = Aj ds J ds' e -(«+«')/* p u (s, s', A), 

where we accepted for the values of the Borel parameters to be the same for the both channels (the issue 
of the heavy quark symmetry): 

ti=t 2 =2t (61) 

and we also suppose that the value of the continuum threshold lo is the same as in two-point sum rules. 
Performing Fourier transformation with respect to A one obtains sum rules for the matrix elements in 
the momentum space: 

^F s 2 tJ D(fc+)e- A/t =4/ ds ds' e -^' 2t Pd (s,s',k + ) , (62) 
z Jo Jo 

-F* t U(l + )e~ A/t =4^ ds J ds' e- (s+s '^ 2t p u (s,s'J + ). (63) 
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The calculation of the diagrams in Fig [2] provides the following analytical results for the spectral densities: 

2 



p u (s,s',l+) = -N c 



1 

47T 2 

1 



1 

Av 



+ 



2 J±_ 

2v 



16 7T 2 Av' 



Av' + \42V, 



5{s) 

2 



+ 

1 m 



[{)<!+< 2s'v' + ] 



1 - — - 



o 

16 At 2 
e~< / 64t2 8{s) 



i±- 9(0 <l + < 2s'v') + s 2 5(s') e-™ 2 ° ^ S ( l -±- - ^ 



16i 



21 



1 m u 
16 it 2 " 



<5(s') <5 



16i 



(64) 



N 



N 
4 



j_y i 

Z^ 2 / 4 

2 

\UU) ' 



9[0<k + < 2s] 



4jV, 
ra) 1 



<5(s') 



i _ .ll^ 

16 At 2 



w 1 + 16i 



16 7T 2 4 



i _ _Ln^ 

16 4t 2 



-|t^(0 < k+ < 2s) + s' 2 5{s)e- m ° / 64t2 S ( k, 



(65) 



16i 



The quantity mo is known as vacuum correlation length and defined as m 2 , = (qg(oG)q) / (qq) ~ 0.8GeV 2 
The diagrams with quark condensate in Fig[5] have been computed using the technique of the non-local 
condensate 22,23. In such approach one introduces vacuum expectation value of the non-local operator 



(0\q(x)[x,0]q(0)\0) ~ (0|gg|0) / dv j(v) e 1 ^/ 4 



(66) 



which has to be understood as a model for the partial resummation of the OPE to all orders. Such treat- 
ment allows to escape the singular 6— function terms which appear in the OPE with the local condensates 
[2"3"j . This is a very general situation and it arises also in the sum rules of the B-meson LCDA [TSJ [TB] . 
For the spectral function we accept the simplest model suggested in [S5J 123] : 



f{v) = 8{v - ml/A), 



(67) 



Let us also remark that we neglect the terms with the gluon condensates because corresponding contri- 
butions are small. The similar observation was made also in the sum rules for B-meson LCDA [TBI [16] ■ 
In the numerical calculations of sum rules (|6"2"]) and (|63p we substitute the value of decay constant F st 
obtained from the corresponding the two-point sum rule [191 120] : 



1 



-F 2 (/i)e- A /* = ^ 
° VW 2tt 2 



ds s 2 e~ 



s/t _ 1 



1 



16i 2 



(68) 



It is instructive to consider the so-called local duality limit t — > oo. Then the sum rules expressions 
are simplified and one obtains 



1 u; 2 



k+\ 



v' + 8tt 2 2lu v' + 



1 - 



2loqv' + J 



- ) 9 [ < 1+ < 2u v' + ] 



(69) 
(70) 



As one can see, the both functions are localized in the region k + < 2lu v + . We expect that this is 
valid only for the leading order approximation, similar to the B— meson LCDA [IB]. Another important 
property is the "good" behavior in the limit k + — > 0. Such behavior at small values of the momentum 
fraction do not contradict to the existence of the convolution integrals (f33|) . 

In the numerical estimates we use for the Borel mass t and continuum threshold ujq the same values 
as in the the two-point sum rules |19j 



0.3 GeV < t < 0.6 GeV, uj q = 0.8 - 1.0 GeV, 



(71) 
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and (uu) — — (240MeV) 3 , m 2 , ~ 0.8GeV 2 . From the expressions for the spectral densities (|64p and (|65|) 
one can easily find that 

pnf = J dk+jya^ U nf = J |±[/nf (fc+) = _xpf /t / +j (72 ) 

with 1/ = Mb/Mjj = 2.63. Therefore we provide the results only for one quantity V nl . Numerically we 
obtained 

V nl = (1.32 ± 0.16) x 10" 2 GeV 2 , (73) 
where the uncertainty arises from the t— and wo— variations. Hence assuming (|72|) we obtain 

a ni „ CiV nf I (g . . ^ ^ - ^ + ie x<jp e*"e^ [e d + ^) } ■ (74) 

Substituting the leading order value for the coefficient function C\{vnh = 4.8GeV) = 1.12 one has 

10 3 a nf ~ (e; • e* D ) (-8.64GeV 2 ) + iex ap s*°e% (-1.20GeV 2 ) . (75) 

Comparing this result with the analogous expression for a 1 (|53j) we observe that both form factors a f 
and a , from the factorizable and non-factorizable soft matrix elements are of the same order. From the 
structure of expressions (|53|) and (f75j) it is easy to see that factorizable contribution dominates in the 
physical form factor F\ ~ 3.8 x 10 _9 GeV but the non-factorizable term provides the largest contribution 
to the the second physical form factor F2 ~ —2.5 x 10 _9 GeV. 

3.3 Branching fraction estimate and conclusions 

With the above results we can estimate the branching fraction. Using for the CKM matrix elements 
I Kid I = 0.974 and \V c b\ = 0.415 and a = 1/137 we obtain for the branching ratio 

B(B° -> D*°j) = (1.52 ±0.35) x 10~ 7 (76) 

The uncertainty given in (|76|) originate from the uncertainties of the hadronic matrix elements. We 
observe that our estimate is of two order magnitude smaller that the experimental bound B(B° — > 
D*°j) < 2.5 x 10~ 5 9 . Our estimate is also significantly smaller than the values provided by previous 
considerations [101 QJJ [12] . The main conclusion is that such small quantity most probably can not be 
measured at existing £?-factories. 

On the other hand our estimate has to be considered carefully. We have used only the leading order 
contribution. There are a lot of corrections which a priory may be of considerable size. We did not 
consider the resummation of the possible large Sudakov logarithms associated with the choice of the 
factorization scale. Typically, the large corrections arise also from the next-to-leading contributions to 
the jet functions in the SCET approach because the corresponding hard scale ~ KE 1 is not very large 
|24j . The complication also arises due to the fact that there are two different heavy quark masses mj 
and m c that introduce an additional scale ambiguity. Therefore on the background of these remarks 
our result has to be considered only as a leading order qualitative estimate. However, we expect that 
all effects mentioned above can not provide such strong enhancement that can make the value of the 
branching measurable for BABAR or BELLE experiments. 
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